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ABSTRACT 

Using the IRAC images from the Spitzer c2d program, we have made a survey 
of mid-infrared outflows in the p Ophiuchi molecular cloud. Extended objects 
that have prominent emission in IRAC channel 2 (4.5 /an) compared to IRAC 
channel 1 (3.6 pva) and stand out as green objects in the three-color images 
(3.6 /iin in blue, 4.5 pm in green, 8.0 pm in red) are identified as mid-infrared 
outflows. As a result, we detected 13 new outflows in the p Ophiuchi molecular 
cloud that have not been previously observed in optical or near-infrared. In 
addition, at the positions of previously observed HH objects or near-infrared 
emission, we detected 31 mid-infrared outflows, among which seven correspond 
to previously observed HH objects and 30 to near-infrared emission. Most of the 
mid-infrared outflows detected in the p Ophiuchi cloud are concentrated in the 
L1688 dense core region. In combination with the survey results for Young Stellar 
Objects (YSOs) and millimeter and sub-millimeter sources, the distribution of 
mid-infrared outflows in the p Ophiuchi molecular complex hints a propagation 
of star formation in the cloud in the direction from the northwest to the southeast 
as suggested by previous studies of the region. 

Subject headings: infrared: ISM — infrared: stars — ISM: individual (p Ophi- 
uchi) — ISM: jets and outflows — stars: formation 



1. Introduction 



Mass outflows play an essential role in the proc ess of star format i on and have been found 
to be ubiquitous in various stages of star formation (lArce et al.ll2007l ; iBally et al.ll2007l ). The 
specific angular momentum of a star-forming mo lecular core is abou t four orders of magnitude 
times higher than that of a typical T Tauri star (jBodenheimerlll995l ) . One way to transfer the 
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excessive angular momentum from the star-forming cores so that the circumstellar m aterial 
can be accreted onto the central stars is through mass outflows ( IShang et al.l 120071 ). Mass 
outflows from YSOs have been observed at different wavelengths. In the visual, high velocity 
jets and HH objects with a typical veloci ty of 100 - 300 km s" 1 trace material ejected by 
the star or shocked surrounding medium (IReipurth fc Ballyll200ll ). Near- infrared molecular 
hydrogen emission features, which usually have veloc ities of several tens k m s , trace H2 gas 
in the jets or in the entrained surrounding medium (jEisloffel et al.ll2000l ). In the millimeter 
wavelength, hig h velocity CO g as with velocity in the range from a few to ten km s" 1 have 
been observed (IWu et al.l I2004J ) . The high velocity CO gas probes the swept or entrained 
medium. The different velocity ranges and excitation conditions for HH objects, near-infrared 
molecular hydrogen emission, and CO outflows indicate the transfer of energy and momentum 
from the high velocity material to the surrounding medium. 

Shocked gas in outflows has abundant atomic and molecular hydrog en line emission 
i n the mid-infrared. Hydrodynamic simulations of outflows from YSOs by [Smith fc Rosen 
( 120051 ) have shown that outflows have strong H2 line emission in all the IRAC four bands. 
In fact, mid-in frared observation s with IRAC aboard of Spitzer have discovered a lot of 
new outflows (TMerm et ah 2008h and plenty new features of prev iously known outflows 



( INoriega-Crespo et al. 



2004 



Smith et al.l 120061 ; iTeixeira et al.l 120081 ) , showing that Spitzer 



IRAC imaging is a powerful tool in the survey of outflows, particular for those that are 
deeply embedded in clouds due to the much less extinction in the mid-infrared compared 
to the optical and near-infrared wavelengths. The Spitzer legacy program c2d performed a 
complete imaging of five nearby star-forming regions with the IRAC and MIPS instruments 
( lEvans et al.ll2003l . 120091 ). which provides a valuable data set to search for mid-infrared out- 
flows in these regions. In this paper we present the results of a mid-infrared outflow survey 
of the p Ophiuchi molecular cloud based on the c2d archive images. 

The p Op hiuchi molecular clo ud is one of the closest star-forming regions at a distance of 
about 130 pc (jWilking et al.ll2008l ). It is intermediate in star formation activity compared to 
the isolated star formation in the Taurus cloud and the rich clusters in the Orion cloud. The 
large scale str ucture of the p Ophiuchi cloud has been revealed by extensive molecular 13 CO 
line mapping ( Lorenlll989al ibl). Millimeter and submillimeter surveys of the central part of 



the cloud have revealed the clumpy structure of t he cloud and have f ound about 55 starless 
cores with masses in the range of 0.02 to 6.3 M (IMotte et al.lll998l ; I Johnstone et al.ll2000l ; 
Young et al.l 120061). The cloud has been surveyed for young stars a t wavele ngths from X 



rav rtMontmerle et al.l 11983c Imanishi et al 



ter (jStanke et al 



198?{~ 2005h . near-infrared (Greene fc Young 



2001 



2006 



1992 



Ozawa et al. 



Al 



(TYoung et al. 1986 ; Wilking et al. 2001 ; Padgett et al. 20081 ). to millimeter a nd submillime- 



en et al. 



2005), visual (jWilking et al. 



20021 ). mid- and far- infrared 



Andrews fe Williams! 120071 ). Based on these surveys, IWilking et al. 
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( 120081 ) compiled a list of 316 young stars in 



About 46 HH obje cts (IWilking et al. 



emission features (IGrosso et al. 



velocity CO outflows (IBontemps 



2001 



1997 



1996 



1688, the main cloud of the p Ophiuchi complex 



Wu et a 



Gomez et al. 



2002 



2003 



Bussmann et al 



Phelps fc Barsonyl 120 04). 119 H 2 



Khanzadvan et al 



200J), 16 high 



20071 ) have been observed in the p 



O phiuchi complex. For a summary of the p Ophiuchi complex we refer to the recent review 
bv lWilking et all (hoosh 



2. The Data 



2.1. The Spitzer Data 

The IRAC observations of p Ophiuchi were conducted on 2004 March 28-30 and 2005 
September 15 and 18. These observations are part of the Spitzer legacy program c2d 
( lEvans et al.ll2003l ) and the program ID is 177. Fig. [1] shows the three-color image of the area 
mapped with IRAC which is about 8.0 d eg 2 . This area is included in the region mapped with 
MIPS (see Fig.l of iPadgett et al.l 120081 ). Each region in Fig. [T] was individually mapped in 
two epochs with an exposure time of 12 seconds (24 seconds in total) and one map was made 
i n the high dynam ic range (HDR) mode, which added a short exposure time of 0.6 seconds 
( lEvans et al.ll2003l ). Six off-cloud fields were observed in order to sample the background 
source counts. 

We used the final mosaic images and source catalogs of the fourth delivery of IRAC and 
MIPS data. The data were processed with Spitzer Science Center pipeline S13 and calibrated 
by c2d team (details can be found in the Delivery Document^). In the final mosaic images 
bright source artifacts have been removed. However, some instrumental effects still survive. 
For example, optical ghosts produced by internal reflections within the filters (labelled with 
"gh" in the subsequent figures) and stray lights produced by illuminating sources off the 
edges of the detector arrays (labelled in the subsequent figures with "Stray light") remain 
in the final mosaic images (see the IRAC data handbook^ for details). 

A catalog of young stellar object candidates in the p Ophiuchi, YSOc, is also delivered 
by the c2d team. lEvans et al.l (120091 ) improved over this catalog. They removed some sus- 
pect sources, added known sources and data obtained at other wavelengths, and calculated 
additional quantities such as bolometric temperature and bolometric luminosity. Their list 



1 http:/ /data. spitzer. caltech.edu/popular/c2d/20071101_enhanced_vl/Documents/ 
c2d_del_document.pdf 

2 http:/ /ssc. spitzer. caltech.edu/irac/dh/iracdatahandbook3. 0.pdf 
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makes a relatively pure YSO candidate sample which includes 292 YSOs. We used their 



l ist of YSOs in this paper. Two schemes for YSO classification were used in lEvans et al. 
( 120091 ). One is the spectral index that was obtained using a least-squares fit to all pho- 



tomet ry between 2 /im and 24 pm, and the criteria for classification are f rom I Greene et al. 



(1994). The other sc heme used the calculation of bolo metric temperature (IChen et al.lll995l ; 



Dunham et al.ll2008l ). Using the calculated spectral indices or bolometric temperatures, 292 



YSOs in the cloud can be classified into 4 classes: 27 in Class I or Class 0, 44 in the flat- 
spectrum Class, 179 in Class II, and 42 in Class III. Note especially that there are three 
Class sources in Ophiuchus, IRAS 16293-2422, MMS 126, and VLA 1623-243, which were 
distinguished from Class I sources on the basis of their bolometric temperatures. 



2.2. The FCRAO Data 



To analyze the distribution of detected outflows in the p Ophiuchi molecular complex, 
we used the 13 CO 1-0 (110.201GHz) emission data from the FCRAO 14 m telescope at New 
Salem, MA, U.S. A. The data were obtained during 2002-2005 with the On-the-Fly (OTF) 
mapping method (IRidge et al.ll2006l ) and are part of the COMPLEtM survey of star-forming 
regions. These 13 CO observations cover 10.0 deg 2 of the p Ophiuchi molecular complex and 
the FCRAO Half Power Beam Width (HPBW) measures 44". The 13 CO integrated intensity 
maps that have been re-gridded to 23"were used. In our 13 CO emission contour maps, the 
units of integrated intensity is antenna temperature (T^) in Kelvin. 



2.3. Identification of IRAC Mid-infrared Outflows 



IRAC channels cover a relatively broad wavelength range from 3.6 pm to 8.0 pm. In 
this wavelength range there are many H 2 lines and atomic and ionic fine-structure lines. 
Rosenthal et al.l (120001 ) have shown near- and mid-infrared spectra of the brightest H 2 emis- 
sion peak of the Orion OMC-1 outflow which are taken with ISO-SWS. The spectra revealed 
a wealth of molecular lines of H 2 , CO, and H 2 0, and a number of H I recombination lines, 
atomic and ionic fine-structure lines. The IRAC 4.5 pm b and contains stron g H 2 lines 
( v = - 0, S(9, 10, 11)) and CO bandheads (v = 1 - 0) jReach et al.l l2006h . Like the 
near-infrared H 2 emission at 2.12 pm, the emission of these lines can be exc i ted by the inter- 
action of protostellar outflows with the ambient medium. ISmith fc Rosen! (120051 ) presented 
synthetic Spitzer IRAC band maps from hydrodynamic simulations of processing protostellar 



3 http: / /cfa- www.harvard.edu/COMPLETE 
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jets and they predicted that the integrated H2 line luminosity is greatest in IRAC band 2 
(4.5/im), about an or der of magnitude higher than those in other IRAC bands (see Table 3 in 
Smith fc Ro sen 200 a) . Observat i ons by Spitzer IRAC h ave in general confirmed this predic- 



tion (Noriega-Crespo et al.ll2004l ; Ijprgensen et al.ll2006l ). According to iNoriega-Crespo et al. 



(120041 ). the remaining emission in the Spitzer IRAC band 2 image after the subtraction 
of the emission in band 1 is mainly from shocked H o and CO gas. Us i ng th e three-color 



images constructed from IRAC bands 1, 2, and 4., iCyganowski et al.l (120081 ) compiled a 
catalog of outflow candidates with the Spitzer Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE) data. They identified the Extended Green Objects (EGOs) in 
the three-color images as outflow candidates. In this paper, we use the difference image con- 
structed from IRAC channel 2 (4.5 /mi) minus channel 1 (3.6 /im) to search for mid-infrared 
outflows. We use the three-color (3.6 /mi in blue, 4.5 /im in green, 8.0 /xm in red) images to 
confirm that the identified outflows are EGOs. 

Images covering the common sky area of the p Ophiuchi cloud in the four IRAC bands 
were trimmed from the downloaded c2d archive mosaic images. The trimmed images were 
then rotated so that the image north direction is top and the east is to the left. For the 
subtraction of emission in band 1 from that in band 2, we selected 80 bright stars in the 
images that are roughly evenly distributed over the sky and are not saturated in either band 
1 or band 2. The average flux ratio in band 1 to band 2 of these 80 stars is used to adjust 
the flux level of stars in band 1 to that in band 2. The resultant difference images were 
visually inspected to search for mid-infrared outflows. We selected extended sources that 
have significant excess at 4.5/im as mid-infrared outflows. Tow ard the core of L 1688, several 



surveys for HH objects and Ho emission have been carried out (IWu et al.ll2002l ; iGomez et al. 



20031 ; iKhanzadyan et al.ll2004l ) and 46 HH objects, including components of HH objects, and 
119 H 2 emission features have been reported. We examined whether there are mid-infrared 
counterparts of these HH objects and H 2 emission features at their corresponding positions 
in the IRAC images. 



3. Results 



We discovered 13 IRAC mid-infrared outflows. Their coordinates and figure reference 
are list in Table [TJ The likely exciting sources for these outflows and their Lada classes are 
given in Table [2J Note t hat IRAS 16293-2 422 is classified as a Class source based on its 
bolometric temperature (jEvans et al.ll2009i ). We also identified 31 mid-infrared counterparts 
of previously known HH objects and H 2 2.12 /mi emission features. The coordinates of these 
mid-infrared counterparts of the previously known outflows, together with the corresponding 
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HH objects and H2 2.12 /im emission features, are presented in Table [31 In this paper we 
designate all the detected IRAC mid-infrared outflows as EGOs. 

The EGOs are numbered by order in right ascension. If the overall morphology of a 
group of knots suggests physical association among them, the knots are considered as parts 
of one EGO and the individual knots are distinguished with additional letters to the EGO 
numbers. On the other hand, if the knots in a region are relatively separated from each other 
and the relationship among them is not clear, each knot is treated as an individual EGO. 
Apparent sub-structures in a knot are labelled with additional numerals, such as EGO 21bl. 

The locations of all the mid-infrared outflows identified in the p Ophiuchi cloud are 
shown in Fig. [2j The grey scale image in Fig. [2] is the IRAC channel 2 (4.5 /im) image. 
The outflows newly discovered in this work are marked with circles and the counterparts of 
known flo ws are labelled with p luses. Two YSO aggregates in the region, L1689- and L1709- 
aggregate (IPadgett et al.ll2008l ). are marked with open squares. The identified mid- infrared 
outflows are located mainly in three regions, i.e., the L1688 core, the IRAS 16293-2422 
region, and the L1709 YSO aggregate region. The images of the newly discovered IRAC 
mid-infrared outflows are shown in Figs. [3TTT71 and the images of mid-infrared counterparts 
of HH objects and H2 emission features are presented in Figs. [T8"ll4"0l 



3.1. New IRAC Mid-infrared Outflows 

Thirteen new mid-infrared outflows have been discovered in our survey. For each mid- 
infrared outflow, we present the channel 2 (4.5 /im) minus 1 (3.6 /im) image and the three- 
color image. 

Figures 3-4 show the region of objects EGO 04 and EGO 12. The distance between the 
two objects is about 3'. Each object consists of two knots. EGO 04 and EGO 12 are only 
detected in IRAC channels 1 (3.6 /im) and 2 (4.5 /im). It is brighter in channel 2 (4.5 /im), 
therefore, appears as an EGO in the three-color image (Fig. 0]). There are three YSOs to the 
sout heast of EGO 04 a nd EGO 12, among which we suggest BKLT J162624-241616 (Class 



II in lEvans et al.ll2009i ) as the exciting source of the EGO 04 and EGO 12 outflow, as this 
source is also detected in the MIPS 24 /im image. The distance from BKLT J162624-241616 
to EGO 04 is 4!7 and the distance to EGO 12 is 1(8. 

Outflow EGO 33 (Figs. 5 and 6) is a diffuse nebula and is visible in IRAC channels 
2 (4.5 /im), 3 (5.8 /im), and 4 (8.0 /im). Mid-infrared co unterparts of the known near- 



infrared outflow in this region, three H 2 emission features in iKhanzadyan et al.l (120041 ). are 



also detected in the IRAC images. We designated these mid-infrared counterparts as EGO 
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34a-c (see Table [3]). EGO 33 appears green and yellow in the three-color image (Fig. [6]), 
which indicates that EGO 33 may have so me emission from dust or PAH. There are two 
YSOs in th is region, MMS12 6 (Class in Istanke et~aD l2006h and BELT J162816-243657 
(Class II in lEvans et aDl2009h. BELT J162816-243657 is associ ated with 1.3 mm and mid- 



infrared emission (lAndre fc Montmerleill994! ; iGreene et al.lll994j ). These YSOs are visible in 
all IRAC four channels and in the MIPS 24 fim image. We suggest BELT J162816-243657 
as the driving source of EGO 33 on the basis that MMS126 has close association w ith the 
known outflow detected in near-infrared and in CO (3-2) emission (IStanke et al.ll2006l ). EGO 
33 is unlikely to be associated with MMS 126 because EGO 33 is not coincident with the axis 
of the known CO outflow. The distance from BELT J162816-243657 to EGO 33 is about 
V.9. 

Both EGO 35 and EGO 36 are complexes of emission filaments ( Figs. 7 and 8). Their 
components are list in Table [TJ EGOs 35 and 36 are visible in IRAC channels 1-3 and are 
brightest in chann el 2 (4.5 iim). Th e nearby YSOs are shown in Fig. [9j BELT J162908- 
241549 (Class II in lEvans et alJl2009h is the nearest YSO to EGOs 35 and 36. However, the 
driving source of EGOs 35 and 36 should be located in the direction of the line connecting 
EGOs 35 and 36 if there is some association between EGOs 35 and 36. Other YSOs in this 
region are located to the northeast or southwest of EGOs 35 and 36, all roughly aligned 
with EGOs 35 and 36. The morphology of EGO 35 and EGO 36 is bow-shock-like, with 
wings extending to the southwest, which suggests that the outflow source should be located 
to the sout hwest. HH 677 is a bout 30" to the northeast of SR 10 and it may be driven by 
this source (IGomez et al.l 120031 ). We can see that EGOs 35 and 36 and HH 677 are roughly 
located on a line. Therefore, it is possible that they are physically associated. 



EGOs 37-40 are located in the L1709 YSO aggregate (jPadgett et al.ll2008f ) (Figs. 10- 
11). EGO 37 is a faint knot and EGO 38 is a faint bow structure and they may have some 
connecti on. They are v i sible in IRAC channels 1 and 2. GWAYL 4, a l ow- luminosity class I 



source m 



Greene et al.l (119941 ) and classified as a flat spectrum source in lEvans et al.l (120091 ) 



is located in the northeastern prolongation of the EGOs 37-38 outflow. The distance from 
GWAYL 4 to EGO 37 and EGO 38 is about 3^6 and 2?6, respectively. On the basis of the 
alignment we suggest GWAYL 4 as the exciting source of EGOs 37 and 38. EG Os 39 and 40 



are lo cated to the southeast and east of YSO L1709-3, a flat-spectrum source in lEvans et al. 
J2009h . EGO 39 and EGO 40 are visible in all IRAC four channels and EGO 39 is also 



visible in the MIPS 24 /xm image. The L1709 YSO aggregate is a small cluster of six YSOs 
and it is difficult at present to identify which YSOs in this cluster as the exciting sources of 
EGO 39 and EGO 40. 



EGO 41 consists of a series of knots, EGO 41a-i, and EGO 42 is a relatively isolated 
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patch (Fig. [T2l . Knots EGO 41a-i constitute a bowlike structure facing to northeast. The 
bowlike structure is visible in all IRAC four channels while EGO 42 is only visible in IRAC 
channel 2. The symmetry axis of the bowlike structure suggests that the driving source for 
EGO 41 is located in the NE-SW direction. SSTc2d J163145.81-243909.0, a Class II source 



( lEvans et al.ll2009l ). is located about 8 ' to the southwest of EGO 41 (see Fig. [2]). This source 
is a possible driving source of EGO 41. 

E GO 41 is located about 5'2 to the southwest of the Class object IRAS 16293- 
2422 (lEvans et al.l 120091 ) which drives two mo lecular outflows, one in the east- west direc 



tion and another one in the NE-SW direction (jWootten &: Lorenlll987l ; ICastets et al.l 12001 



Hirano et al.ll200ll ; iLis et al.ll2002l ; IStark et al.ll2004h . The IRAS 16293-2422 E-W molecular 
outflow also has its spectacular demonstration in the mid-infrared which is designated as 
EGO 43 and, together with EGOs 41 and 42, is shown in Fig. [TJ The scale of the EGO 
43 outflow measures about 11'. It consists of four main portions, EGO 43 A-D. The details 
of the inner portions A-B and the outer portions C-D of this outflow are sho wn in Figs. [T41 
andll5l respectively. IRAS 16293-2422 is a binary with components A and B ( Mundy et al.l 
19921 ) and their positions are marked in Fig. [TH with pentagrams. IRAS 16293-2422 is vis- 
ible in the MIPS 24 /xm image and invisi ble in the IRAC im ages. In Fig. dH the position 
of the submillimeter source IRAS 16293E (jNutter et alil2006h . 90" to the southeast of IRAS 
16293-2422, is also marked. 

From Fig. [12] we can see that the symmetry axis of the EGO 41 bowlike structure goes 
through EGO 42 and IRAS 16293-2422. This axis c oincides well w ith the axis of the IRAS 
16293-2422 NE-SW molecular outflow (see Fig. 12 in lLis et alJliooi ). suggesting that EGOs 
41 and 42 may have physical associatio n with this molec ular outflow which is attributed to 
the A component of IRAS 16293-2422 Jstark et al.l 120041 ). 



We note that the EGO 41 bow shock points to source IRAS 16293-2422. This morphol- 
ogy usually suggests that the driving source is located on the direction opposite to IRAS 
16293-2422. However, bow shocks are not always in a c onfiguration of facing away from their 
driving sources. For example, the HH 311 b ow shock (IReipurth et al.l 119971 ) and the south- 
eastern bow shock of the V380 Ori-NE flow (IDavis et al.ll2000l ) are in a configuration facing 
toward their exciting sources. This configuration can be under stood if the flow impacts a 
dense, ambient clump, around which shocked gas is streaming ( IDavis et al.l 120001 ) . So it is 
also possible that EGO 41 is associated with IRAS 16293-2422. 

EGO 44 is a bow shock in the IRAC channel 2 minus channel 1 image (Fig. [TBI) . It is 
visible in all the IRAC channels and Fig. [T7] presents the three-color image of the region. 
EGO 44 is located at the east outskirts of the p Ophiuchi cloud (F ig.[2]). The nearest known 
young stellar object to EGO 44 is IRAS 16367-2356 (Class II in lEvans et al"1l2009h which 
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is about 15' to the west of EGO 44. As the EGO 44 bow shock faces to the east, IRAS 
16367-2356 is likely the exciting source of EGO 44. 



3.2. IRAC Counterparts of Known Outflows 



Extensive surveys of outflows toward the p Ophiuchi cloud, including HH objects in the 
optical and H 2 2.12 /xm emission in the near-infrared, have been performed. In total, 46 HH 
objects, inclu ding components of HH objects, and 119 H 2 emission f eatures have be detected 
in this region ( 



Grosso et al.ll200ll : IWu et al.ll2002l : iGomez et al.ll2003l : iKhanzadvan et al.ll2004 



Phelps fc Barsonyl 120041 ) . We identified in the IRAC images 31 EGOs that correspond to the 
known outflows, among which seven EGOs correspond to known HH objects and 30 EGOs 
to H 2 near-infrared emission. 

Figs. [TBI and [T91 show the region of EGO 01 which has been identified at 2.12 /iin 
by iLucas et al.l (120081 ). EGO 01 are three diffuse nebulae in the 2.12 /xm image while it 
consists of 3 faint knots in the 3.6 jttm and 4.5 pm images. Two YSOs, GSS 23 and SR 
4, are located in the nearby of EGO 01. GSS 2 3 is c lassified as a weak- line T Tauri star 
(IBouvier fc Appenzellerl Il992l ) . iBary et al.l (120021 ) and iBitner et al.l (120081) detected narrow 
H 2 e mission surrounding GSS 23. SR 4 is an emission-line star (jstruve fc Rudkjobing 1949 ) 
and iPhelps &: Barsonvi (120041 ) detected an HH object, HH 312, to the southeast of SR 4. 
Lucas et al. ( 20081 ) suggested that there are some connection between the H 2 emission and 
the outflow driving by VLA 1623-243 (see Figs. [241 - [261 and discussions on EGOs in that 
region) on the basis of their locations and morphology. 

Figs. 5-6 have shown the region of E GO 34 which are counterpa rts of the outflow 



[KGS2004] f05-04 (IKhanzadvan et al.l 120041 ) . IKhanzadvan et all (l2004f ) and ISmith et al. 



( 20051 ) observed three knots which correspond to EGO 34a-c in their near-infrared ob- 
servations. However the Spitzer images show much more details of the outflow than the 
near-infrared images. In Figs. [5] and [6] we can see extended diffuse emission to the northeast 
of EGO 34c and to the southwest of EGO 34a. The central source of the outflow, MMS126, 
wh ich is identified as a low-mass Class object in the millimetre continuum observations 
by IStanke et al.l (120061 ) (Classified as Class in lEvans et al.l 120091) is visible in all IRAC 
bands and also in the MIPS 24 /zm image. IStanke et al.l (120061 ) has observed a molecular 
CO outflow from MMS126 which is orientated in the NE-SW direction. 



Figs. [20] and [21] present the images of EGOs 02, 06, 08, 15, and 16 that correspond s 
to near-infrared fl o ws [K GS2004] f 10-01 and [KGS2004] 10-02 fchanzadvan et al.l l2004h . 
Khanzadyan et al.l (120041 ) detected nine H 2 emission components of [KGS2004] fl0-01. How- 
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ever, only four knots, EGOs 02, 08, 15, and 16, can be detected in the IRAC images. 
These four EGOs a re located on the two side s of BKLT J162619-242820, a Class II object 
favans et al1l2009h . iKhanzadvan et all (120041 ) proposed BKLT J162619-242820 as the driv- 
ing source of this flow. EGO 06 has two components, EGO 06a and 06b, corresponding to 
[KGS2004] flO-02. The driving source of EGO 06 is unclear. 



The IRAC images of EGOs 03, 05, and 10 are presented in Figs. [22] and El EGO 03 
is the counterpart of [GSWC2003] 8a a rid EGO 05 is the c ounterpart of [ GSW C20031 8d. 
EGO 10 corresponds to [GSWC2003] 9 JGomez et alil2003h . iGomez et all hoO$ \ identified 
[GSWC2003] 8a as three di stinct components . In Figs. [22] and [23] EGO 03 is a diffuse patch. 
GSS 26, a Class II source fevans et alil2009h . is the nearest YSO to EGO 03 and EGO 05 
and may be the driving source for these EGOs. EGO 10 is a diffuse structure within which 
we can identify several relatively bright knots. The components of EGO 10 constitute a faint 
bow-like structure with its symmetry axis pointing to GSS 26, suggesting that EGO 10 is 
also driven by this source. 

Fig. M shows the locations of EGOs 07, 09, 11, 13-14, and 17 -18. Their p r evious 
identifications in optical and near-infrared can be found in Table [3] I Gomez et al.l (120031 ) 
observed this region and proposed that most of the emission featu res in this region are 
driven by VLA1623-243, a Class source in L1688 (jAndre et al. 1993 ) which is marked with 
a pentagram in Fig. [2H VLA1623-243 has been identified a s the driving sour ce of a large 
scale m olecular outflow which lies in the direction of NW-SE flAndre et al.lll990l ). iDent et al. 
(119951 ) presented maps of the large-scale 12 CO outflow and H 2 near- infrared emission images 
of this region. 



Fig. [25] shows the northwest part of the VLA1623-243 region. EGOs 07, 11, arid 13 - 
14 distribute roughly as a chain in the NW-SE direction. ICaratti o Garatti et al.l (120061 ) 
observed this region in near-infrared and obtained proper motion measurements for the 
knots in the VLA1623 region. EGO 11 correspon ds to HH 31 3 B, wh ich is likely associated 
with GSS 30, a class I binary system candidate in lChen et al.l (120071). as the positio n angle 
of proper motion measured for HH 313B is about 24 3°(Caratti o Garatti et al.l 120061 ) . EGO 



09 is a diffuse nebula on the bottom right of Fig. [25] IKhanzadvan et al.l (120041 ) detected this 



feature in the near-infrared and suggested that it is driven by a source to the north. 



Fig. [26] shows the VLA 1623-243 southeast region. We can see a chain of knots in 
the EGO 17 and EGO 18 flows. EG O 17 and EGO 18 correspond to near-infrared flows 
[GSWC2003] 20 and 21, respectively (jGomez et al.ll2003l ). The bright mid-infrared emission 
knots are visible at the IRAC 3.6 ^m, 4.5 /iin, and 5.8 /iin images while the faint ones are 
only visible in the 3.6 /iin and 4.5 jttm images. 
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The IRAC images of EGOs 19, 21, and 25 are shown in Figs.[27]and [28l EGO 19 is a bow- 



like structure which corresponds to [GSWC20 03] 24a in the near-infrar ed ([Gomez et al. 



20031) 

and an HH object candidate in the optical (jPhelps fc Barsonyll2004l ). IGomez et al.l (j2003l ) 
suggested Elias 2-26 (Class II in lEvans et ali l2009h as the driving source of [GSWC2003] 24a. 
The distance from Elias 2-26 to EGO 19 is about 2'2. EGO 21a is a faint patch and EGOs 
21bl and 21b2 are elongated knots. EGOs 21bl and 21b2 correspond to HH 79al and a2 in 



the optical image (jPhelps fc Barsonyll2004l ). The relationship between EGO 21a and EGO 



21b is uncertain and the y could be excited by different exciting sources (IGomez et al.l 12003 



Phelps fc Barsonyll2004j ). EGO 25 consists of two bright knots on the top left o f Figs. [271 and 



281 T hev correspond to [GSWC20 03] 25a and 25b in the ne a r-infrared ima g es (I Gomez et aL 



20031 ) and HH 711 in the optical flPhelps fc Barsonvl I2004T ). IGomez et al.l (120031 ) suggested 



three possible driving sou rces for this outflow, incl uding; BKLT J162658-241836, VSSG 24, 
and SR 21. In contrast, iPhelps fc Barsonvl (120041 ) proposed that VSSG 3 is the exciting 
source of HH 711. 



Figs. [29] and [30] show the region of EGOs 20 and 22-23. EGO 20 is a faint knot on 
the bottom right of Figs . [291 and [301 EGO 20 corresponds to a known HH object, HH 314 
JPhelps & Barsonvl 120041 ). EGO 22 is a diffuse knot which c orresponds to [KGS2004] f04-01 
in the near-infrared and HH 673 (jPhelps fc Barsony! 120041 ) in the optical. EGO 23 has a 
bow shock shape and coincides to [ KGS2004] f04-01b. There are three YSOs i n this region, 
BKL T J162641-244015 (Class II in lEvans et al. hood) . BBRCG10 (Class I I inlEvans et al. 
20091 ). and BBRCG19 (Class I, Iwilking et al.l Il989h . IPhelps fc Barsonvl J2004|) proposed 



BKLT J162641-244015 as the driving source of HH 314. 



Figs. I3T1 and [321 show the region of EGO 24, which consists of components EGOs 24a-d. 
These components co rrespond to four emission knots in the near-infrared flow [KGS2004] 
f09-01. In this region iKhanzadyan et al.l (120041 ) detected seven H 2 emission knots, f09-01a- 
g. EGO 24a has a bow shock shape and EGO 24d is a bright knot while EGO 24b-c are 
faint diffuse nebula e. EGO 24a and EGO 24c may belong t o the same bow sh ock structure. 
Gomez et all J2003f ) suggested that BBRCG 24 (Class II in lEvans et aDl2009l ) is the driving 
source of the e mission feature corr esponding to EGO 24d and that WL 16, a early type 
embedded star ([Padgett et al.ll2008l ) located about 3' so uth of EGO 24a, might be responsible 
for the emission features corresponding to EGO 24a-c. IKhanzadyan et al. (120041 ) suggested 
another possibility that f09-01a-g belong to the same flow. 



Figs.[33]and[3lshow the images of EGO 26 that coincides to [KGS2004] f03 -02 fehanzadvan et al 
2004 ) in the near- infrared. EGO 26 also corresponds to an HH object, HH 224S JPhelps fc Barsonvl 
20041 ). BKLT J162722-244807, a Class II source in lEvans et all (I2009I ). is the nearest pro- 
tostar to this EGO. However, considering the association between [KGS2004] f03-02 and 
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other flows in this region, neither iKhanzadyan et al.l (120041 ) nor iPhelps fc Barsonyl (120041 ) 
sug gest this object as the e xciting source of the emission feature corresponding to EGO 
26. IPhelps fc Barsonyl (120041 ) detected several [SII] emission features named HH 224N, HH 
224NW1 and HH 224 NW2 to the northwest of HH 224S. They believed that HH 224S 
belo ngs to a large flow drive n by an embedded source located to the northwest of EGO 
26. IPhelps fc Barsonyl (120041 ) suggested GY 193, a YSO located in the north west about 7 
' away from HH 224S, as the likely exciting source. IKhanzadyan et al.l (120041 ) detected the 
H 2 emission counterpart of HH 224S, but they did not detect the counterparts of HH 224N 
and HH 224NW. Therefore, they suggested EM*SR 24, which is about 8 ' to the north- 
west of [KGS2004] f03-02, as the exciting source based on the fact that [KGS2004] f03-02, 
[KGS2004] f03-01, and EM*SR 24 are aligned along the northwest-southeast direction. In 
the IRAC images, only the counterpart of HH 224S, i.e., EGO 26, is detected. As EGO 26a 
and EGO 26b form a wide-cavity lobe which faces away from BKLT J162722-244807, we 
propose BKLT J162722-244807 to be the exciting source of EGO 26. Figs. [35] and [36] show 
the reg ion of EGOs 27-29 and E GO 31 is shown in Figs. [37] a nd EM EGOs 27-29 a ll are faint 
knots. IKhanzadyan et al.l (120041 ) proposed YLW 15 (Class I. IWilking et al.lll989l ) to be the 
likely driving source of the near-infrared emission feature corresponding to EGO 29. There 
are abundant YSOs in this region and the nearby YSOs are marked in Figs. [33 and [37] 



Figs. [H and M show the region of EGO 30 and EGO 32. IKhanzadyan et al.l (l2004h 
detected eight H 2 emission knots in this region, f08-01a-h. Some of the non-detected near- 
infrared emission knots are located at the position of the stray light marked in Fig. [3H] or 
near the bright star YLW52. The non-detection of these knots may be due to these effects. 
EGO 30, corresponding to [KGS2004] f08-01a, is a bright arc -like structure . EGO 32a-b are 
two faint knots that coincide to [KGS2004] f08-01g and h. iGrosso et al.l (120011 ) proposed 
that EGO 30 arises from the interaction of two different jets emanating from two Class I 
protostars, YLW52 (see Fig. [3HD and YLW15, which is about 10' to the southwest of EGO 
30. 



3.3. Discussion and Conclusions 



The distribution of mid-infrared out f lows in the p Ophiuchi molecular cloud, tog ether 
with that of HH object s JWu et alibooi bomez et al1l2003l : IPhelps fc Barsonyl l2004h and 
molecular CO outflows (jWu et al.ll2004l ). is shown in Fig. UU The grey scale image in Fig. HTj 
is the IRAC 4.5 /im emission and the contours give 13 CO 1-0 emission in the region. It 
can be seen that mid-infrared outflows, HH objects, and molecular CO outflows roughly 
have the same spatial distribution and they are in general concentrated in regions with high 
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gas density, particularly in the region of the L1688 dense core. However, difference in the 
distribution of mid-infrared outflows, HH objects, and molecular CO outflows are apparent 
in Fig. HU First, a substan tial number of HH objects, 17 out of 46 that are detected in the p 
Ophiuchi molecular cloud (IGomez et al.lll998t IWu et al.ll2002l ; iPhelps fc Barsonyl 120041 ). are 
located in regions of 13 CO 1-0 emission less than 2K, while none of mid-infrared outflows 
and molecular CO outflows are detected in these regions of relatively low gas density. This 
result may indicate that the excitation of mid-infrared and CO molecular outflows require 
more dense surrounding medium than HH objects. Secondly, we can see that in the core of 
the L1688 cloud mid-infrared outflows are distributed in the northern and southern parts of 
the core while there are several CO molecular outflows detected in the middle of the core. 
Using outflows as a tracer of current star formation, we can see from Fig. UT] that the core 
of the L1688 cloud is much more active in current star formation than the core of the L1689 
cloud. 

The distribution of YSO candidates in the p Ophiuchi molecular cloud is shown in 
Fig. H2J YSOs in this region are mainly concentrated in the L1688 cloud core. We note 
that abundant YSOs are located in the region to the northwest of the L1688 cloud core 
though there is little dense gas in this region. In fact this part of the p Ophiuchi molecular 
cloud hosts four Class I and two flat-spectrum sources. However, no mid-infrared outflows 
are detected in this region (see Fig. |4"2]) . The non-detection of mid-infrared outflows in this 
region may be due to that little dense gas remains in this region. 

Several extensive surveys have been made toward the p Ophiuchi molecular cloud for 
millimeter sources (MMSs) and sub-millimeter sources fSMMsl (|Woottenlll989l: IWilson et al.l 



Stanke et al.ll2006uYoung et al 



1999 ; IVisser et ali2002l ; Uohnstone et al.ll2004j ; iNutter et al.ll2006l : 

20061 ). The MMSs and SMMs detected in these surveys are shown in Fig. H3j It can be seen 
that the detected MMSs and SMMs are exclusively distributed in the L1688 cloud core and 
the L1689 cloud core. The L1689 cloud core, hosting about 33 MMSs and SMMs, is compa- 
rable to the L1688 cloud core, which hosts about 31 MMSs and SMMs. However, we have 
seen that the L1688 cloud core hosts much more outflows and YSOs than the L1689 cloud 
core, which may mean that the MMSs and SMMs in the L1689 cloud core are in general 
in earlier evolutionary stages than those in the L1688 cloud cores and that they are mainly 
prestellar cores. The distributions of ouflows, YSOs, and MMSs and SMMs suggest a global 
star formation history in the p Ophiuchi molecular cloud which may be that the star for- 
mation in this region took place first in the northwestern part of the p Ophiuchi molecular 
cloud and most dense gas in this region has been dispersed. The L1688 cloud core is now 
undergoing active star formation while most of the dense cores in the LI 689 cloud are still 
at the prestellar stages. This result is consistent with the suggestion that star formation in 
Ophiuchus is triggered by ionisation fronts and winds from the Upper Scorpius OB associ- 
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ation which is located to the west of the p Ophiuchi cloud (ILoren fc Woottenl Il986l ; lLoren 
1989al : iNutter et alil2006L 



In summary, we detected 13 new mid-infrared outflows in the p Ophiuchi molecular cloud 
that have not been detected previously at other wavelengths. In addition, 31 mid-infrared 
outflows which correspond to previously detected HH objects or near-infrared emission are 
detected as well. Seven of these mid-infrared outflow features correspond to previously 
detected HH objects and 30 to near-infrared emission. Most of the detected mid-infrared 
outflows are concentrated in the dense core of the L1688 cloud, with only twelve of which 
locating outside of this core. In combination with the distribution of YSOs and MMSs and 
SMMs in the p Ophiuchi molecular cloud, the distribution of mid-infrared outflows hints a 
propagation of star formation in the p Ophiuchi molecular complex in the direction from 
northwest to southeast. 
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Table 1. IRAC Mid-infrared Outflows in Ophiuchus 
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16:32:28.6 


-24:28:31 


13,15 


EGO 


43b4 


16:32:29.2 


-24:28:16 


13,15 
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Table 1 — Continued 







a 


8 


Figure 


Object 


(J2000.0) 


(J2000.0) 


Reference 


EGO 


43b5 


16:32:30.1 


-24:28:28 


13,15 


EGO 


43b6 


16:32:33.8 


-24:29:04 


13,15 


EGO 


43b7 


16:32:34.7 


-24:28:33 


13,15 


EGO 


43b8 


16:32:35.4 


-24:29:04 


13,15 


EGO 


43b9 


16:32:36.1 


-24:28:28 


13,15 


EGO 


43bl0 


16:32:36.5 


-24:28:32 


13,15 


EGO 


43bll 


16:32:36.6 


-24:28:56 


13,15 


EGO 


43cl 


16:32:44.7 


-24:29:28 


13,14 


EGO 


43c2 


16:32:45.7 


-24:29:09 


13,14 


EGO 


43c3 


16:32:47.1 


-24:29:09 


13,14 


EGO 


43c4 


16:32:47.2 


-24:28:49 


13,14 


EGO 


43c5 


16:32:48.9 


-24:29:13 


13,14 


EGO 


43dl 


16:32:55.9 


-24:28:13 


13,14 


EGO 


43d2 


16:32:57.5 


-24:28:24 


13,14 


EGO 


43d3 


16:32:58.4 


-24:28:39 


13,14 


EGO 


43d4 


16:33:01.3 


-24:28:54 


13,14 


EGO 


43d5 


16:33:02.9 


-24:28:53 


13,14 


EGO 


43d6 


16:33:04.2 


-24:28:40 


13,14 


EGO 


43d7 


16:33:07.2 


-24:28:27 


13,14 


EGO 


44 


16:40:50.6 


-24:01:21 


16,17 



Note. — Units of right ascension are hours, min- 
utes, and seconds, and units of declination are degrees, 
arcminutes, and arcseconds. 
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Table 2. Likely Exciting Sources of the Newly Detected Outflows 



Distance c 



Object 


Exciting Source a 


Lada Class b 


(arcmin) 


EGO 04 


BKLT J162624-241616 


11 


4.74 


EGO 12 


BKLT J162624-241616 


11 


1.78 


EGO 33 


BKLT J162816-243657 


II 


1.88 


EGO 35 


SR. 10, SR 20W, BKLT J162908-241549, SR 20, BKLT J162923-241359 


II, II, II, III, III 


2.84, 16.94 


EGO 36 


SR 10, SR 20W, BKLT J162908-241549, SR 20, BKLT J162923-241359 


II, II, II, III, III 


1.70, 20.58 


EGO 37 


GWAYL 4 


Flat 


3.6 


EGO 38 


GWAYL 4 


Flat 


2.55 


EGO 39 


L1709 3 


Flat 


0.60 


EGO 40 


L1709 3 


Flat 


0.70 


EGO 41 


SSTc2d J163145.81-243909.0, IRAS 16293-2422 


II, 


8.00, 5.50 


EGO 42 


SSTc2d J163145.81-243909.0, IRAS 16293-2422 


II, 


10.36, 3.00 


EGO 43 


IRAS 16293-2422 





11.35 d 


EGO 44 


IRAS 16367-2356 


11 


14.90 



"•S ource designati o ns us ed in the table correspond to the following opti cal or infrare d surv eys: BKLT is 
from iBarsonv et al.l [U997T) , GWAYL from IGreene et all l ll994h , L1709 from IPadgett et"aU 11200811 and SR from 
IStruve fc Rudkiobind l|l949T l 

b The Lada Clas s is classified ba s ed on the spectral index and bolometric temperature. The criteria are from 
IGreene et ail l|l994l l and lChen et~ZI l|l995h 

c The distances to the nearest exciting source and to the farthest exciting source if there are more than one possible 
exciting source. 

d This is the scale of the outflow. 
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Table 3. Mid-infrared Counterparts of Known HH Objects and H 2 Flows in p Ophiuchi 







a 


8 


Previous 


Figure 


Counterpart 


(J2000.0) 


(J2000.0) 


Identifications 3 - 


Reference 


EGO 


Ola 


16:25:49.8 


-24:22:12 


A b 


18,19 


EGO 


01b 


16:25:50.4 


-24:22:22 


A b 


18,19 


EGO 


01c 


16:25:51.3 


-24:22:16 


A D 


18,19 


EGO 


02 


16:26:03.1 


-24:29:33 


[KGS2004] flO-Olb 


20,21 


EGO 


03 


16:26:06.5 


-24:21:35 


[GSWC2003] 8a 


22,23 




05 


16:26:11.1 


-24:20:33 


[GSWC2003] 8d 


22,23 


EGO 


06a 


16:26:12.7 


-24:31:13 


[KGS2004] fl0-02a 


20,21 


EGO 


06b 


16:26:13.5 


-24:31:10 


[KGS2004] fl0-02b 


20,21 


EGO 


07a 


16:26:15.9 


-24:22:41 


[GSWC2003] 11a 


24,25 


EGO 


07b 


16:26:16.4 


-24:22:45 


[GSWC2003] lib 


24,25 


EGO 


08 


16:26:16.6 


-24:28:35 


[KGS2004] flO-Ole 


20,21 


EGO 


09 


16:26:16.6 


-24:25:39 


[KGS2004] fl0-03 


24,25 


EGO 


10 


16:26:17.3 


-24:18:05 


[GSWC2003] 9 


22,23 


EGO 


11 


16:26:17.5 


-24:23:13 


HH 313B C 


24,25 


EGO 


13 


16:26:19.0 


-24:23:05 


HH 313A, VLA 1623A, [DMW95] H5 


24,25 


EGO 


14a 


16:26:21.4 


-24:23:31 


[GSWC2003] 14f 


24,25 


EGO 


14b 


16:26:21.8 


-24:23:37 


[GSWC2003] 14e 


24,25 


EGO 


14c 


16:26:21.9 


-24:23:26 


[GSWC2003] 14g 


24,25 


EGO 


15 


16:26:25.3 


-24:27:34 


[KGS2004] flO-Olh 


20,21 


EGO 


16 


16:26:28.4 


-24:27:05 


[KGS2004] flO-Oli 


20,21 


EGO 


17a 


16:26:31.1 


-24:25:14 


[GSWC2003] 20a, [KGS2004] fl0-04d, [DMW95] H3 


24,26 


EGO 


17b 


16:26:31.7 


-24:25:17 


[GSWC2003] 20b, [KGS2004] fl0-04e 


24,26 


EGO 


17c 


16:26:32.3 


-24:25:19 


[GSWC2003] 20c, [KGS2004] fl0-04f 


24,26 




17d 


16:26:32.9 


-24:25:22 


[GSWC2003] 20d, [KGS2004] fl0-04g, [DMW95] H2 


24,26 




17c 


16:26:33.7 


-24:25:29 


[GSWC2003] 20f, [KGS2004] fl0-04h, [DMW95] HI 


24,26 


EGO 


18a 


16:26:35.7 


-24:25:27 


[GSWC2003] 21a, [KGS2004] fl0-04i 


24,26 




lob 


16:26:36.6 


-24:25:39 


[GSWC2003] 21b, [KGS2004] fl0-04j 


24,26 


EGO 


18c 


16:26:37.5 


-24:25:40 


[GSWC2003] 21c, [KGS2004] fl0-04k 


24,26 


EGO 


18d 


16:26:38.3 


-24:25:46 


[GSWC2003] 21d 


24,26 


EGO 


18c 


16:26:39.2 


-24:25:54 


[KGS2004] 110-04m 


24,26 


EGO 


18f 


16:26:40.2 


-24:26:07 


new d 


24,26 


EGO 


19 


16:26:38.2 


-24:18:29 


[GSWC2003] 24a 


27,28 


EGO 


20 


16:26:39.0 


-24:40:53 


HH 314 


29,30 


EGO 


21a 


16:26:44.4 


-24:20:07 


HH 79b c 


27,28 


EGO 


21bl 


16:26:46.4 


-24:20:15 


HH 79al c 


27,28 


EGO 


21b2 


16:26:47.2 


-24:20:12 


HH 79a2 c 


27,28 


EGO 


22 


16:26:47.5 


-24:38:16 


HH 673, [KGS2004] f04-01a 


29,30 


EGO 


23 


16:26:59.0 


-24:37:03 


[KGS2004] f04-01b 


29,30 


EGO 


24a 


16:27:08.1 


-24:33:05 


[KGS2004] f09-01a, [GSWC2003] 3f 


31,32 


EGO 


24b 


16:27:08.4 


-24:33:32 


[KGS2004] f09-01d, [GSWC2003] 3d 


31,32 


EGO 


24c 


16:27:08.6 


-24:33:17 


[KGS2004] f09-01c, 


31,32 


EGO 


24d 


16:27:09.5 


-24:33:41 


[KGS2004] f09-01c, [GSWC2003] 3b 


31,32 


EGO 


25a 


16:27:10.3 


-24:17:45 


[GSWC2003] 25a, HH 711 


27,28 


EGO 


25b 


16:27:16.4 


-24:43:56 


[GSWC2003] 25b 


27,28 


EGO 


26a 


16:27:22.3 


-24:49:14 


[KGS2004] f03-02a 


33,34 
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Table 3 — Continued 



a S Previous Figure 



Counterpart 


(J2000.0) 


(J2000.0) 




[identifications 51 


Reference 


EGO 26b 


16:27:22.8 


-24:48:53 


[KGS2004] 


f03-02b, HH 224S 


33,34 


EGO 27 


16:27:28.2 


-24:37:56 


[KGS2004] 


f04-06 


35,36 


EGO 28 


16:27:39.6 


-24:41:50 


[KGS2004] 


f05-01 


35,36,37,38 


EGO 29a 


16:27:39.8 


-24:38:09 


[KGS2004] 


f04-03e, [GSWC2003] 7a 


35,36 


EGO 29b 


16:27:40.2 


-24:37:60 


[KGS2004] 


f04-03f, [GSWC2003] 7b 


35,36 


EGO 30 


16:27:43.2 


-24:31:51 


[KGS2004] 


f08-01a 


39,40 


EGO 31 


16:27:50.4 


-24:40:08 


[KGS2004] 


f05-03 


37,38 


EGO 32a 


16:28:03.9 


-24:30:54 


[KGS2004] 


f08-01g 


39,40 


EGO 32b 


16:28:05.2 


-24:30:40 


[KGS2004] 


f08-01h 


39,40 


EGO 34a 


16:28:17.7 


-24:38:09 


[KGS2004] 


f05-04a 


5,6 


EGO 34b 


16:28:19.4 


-24:37:09 


[KGS2004] 


f05-04b 


5.6 


EGO 34c 


16:28:23.1 


-24:36:10 


[KGS2004] 


f05-04c 


5.6 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination 
are degrees, arcminutes, and arcseconds. 

a Source designations us ed in the tab l e corr espond to the fol lowing optical or infrared sur- 
veys: [GSWC2 0031 is fromlGomez et a.1.1 d2003l l. [KGS2004] from lKhanzadvan et allj2004h . and 
[DMW95] from lDent et alj <1995h ~ 

b This identification is from lLucas et aL (2008) and it is not included in Simbad 

c Also detected in the near-infrared and the identifications are designated with HH object 
numbers. 

d This source is a new feature which has no counterpart of optical or near-infrared flow. 




Fig. 1. — Three-color image of the p Ophiuchi molecular clouds with IRAC band 1 (3.6/iin, 
blue), IRAC band 2 (4.5At m, green) an d IRAC band 4 (8.0/im, red). We marked the main 
Lynds cloud in the region (jLyndall962l ). 
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Fig. 2. — The portion of the p Ophiuchi cloud which includes all the IRAC mid-infrared 
outflows. The grey scale image is the IRAC channel 2 (4.5//m) image. The newly detected 
outflows are marked with circles and the counterparts of previously known HH objects and H 2 
flows are labelled wi th pluses. The open squares represent the L1689 and LI 709 aggregates 
( iPadgett et al.ll2008l ). We also labelled three YSOs with fille d squares for the convenience of 
subsequent discussion (see Section 3.1). SR sources are from lStruve fc Rudkjobingi (Il949i ). 
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Fig. 3. — The region of EGO 04 and EGO 12. The image is the difference i mage of IRAC 
chann el 2 (4.5 fim) minus channel 1 (3.6 /xm). The BKLT sources are from iBarsony et al. 
(|l997h . 
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Fig. 4. — Three-color image of the EGO 04 and EGO 12 region. The image uses IRAC 
channel 1 (3.6 /zm, blue), channel 2 (4.5 /xm, green) and channel 4 (8.0 //m, red). The EGOs 
identified in the region are marked. 
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Fig. 5. — The difference image of the EGO 33 and EGO 34 region. Filt er ghosts are marked 
with "gh". MMS126 (Class 0. Istanke et~al1l2006h and BKLT sources faarsonv et al.lll997h 
are marked. 
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Fig. 6. — Three-color image of the EGO 33 and EGO 34 region. Others are the same as in 
Fig.H 
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Fig. 7. — Region of EGOs 35 and 36. Others are the same as in Fig. [3] 
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Fig. 8. — Three-color image of the region of EGOs 35 and 36. Others are the same as in 
Fig. HI Note that the bow shock morphology of EGOs 35 and 36 is outlined with the white 
lines. 




Fig. 9- 



A large view of the EGOs 35 and 36 region. YSOs in this region are marked. 
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Fig. 10.— The region o f EGOs 37-40. The members of the L1709 YSO aggregate, L1709 
1-6 ( IPadgett et al.ll2008l ). are labelled. Filter ghosts and stray lights are marked. The bow 
shock morphology of EGO 38 is outlined with a white line. Others are the same as in Fig. [3] 




Fig. 11. — Three-color image of the region of EGOs 37-40. Others are the same as in Fig. HI 




Fig. 12. — The region of EGOs 41 and 42. The bow shock morphology of EGO 41 is outlined 
with a white line. Others are the same as in Fig. [3j 
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Fig. 13. — Three-color image of the IRAS 16293-2 422 region. EGOs in the region, 41-42 and 
43A-D, are marked. The GWAYL source is from lGreene et al.l (I1994J ). Others are the same 
as in Fig. HI 
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Fig. 14. — Image of the inner portion of the EGO 43 outflow. The locations of IRAS 
16293-2422 A, B and IRAS 16293E are marked with pentagrams. Others are the same as in 
Fig.S 
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Fig. 15. — Image of the outer portion of the EGO 43 outflow. Others are the same as in 
Fig. El 
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Fig. 16. — The region of EGO 44. The bow shock morphology of EGO 44 is outlined with 
a white line. Others are the same as in Fig. [31 
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Fig. 17. — Three-color image of the EGO 44 region. Others are the same as in Fig. HI 




Fig. 18. — The reg ion of EGO 01, a coun terpart of H2 emission (ILucas et al.l 120081). The 
GSS s ource is from iGrasdalen et al.l (119731 ) and the SR source is from IStruve fc Rudkjobing 
(119491 ) . Others are the same as in Fig. [31 




Fig. 19. — The same as Fig. HJ but for the EGO 01 region. 
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Fig. 20.— Th e same as Fig. [1 but for the EGOs 02, 06, 08, 15, and 16 region. The BKLT 
source is from iBarsony et al.l (119971 ). 
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Fig. 21.— The same as Fig. SJ but for the EGOs 02, 06, 08, 15, and 16 region. 
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Fig. 2 2. — T he same as Fi g. [31 but for the r egion of EGOs 03, 05, and 10. Elias 2 is from 
Eliad (119781 ). GSS 26 from iGrasdalen et al.l (119731 ). and ISO-oph 18 from iBontemps et al. 
(120011 ). The bow shock morphology of EGO 10 is outlined with a white line. 
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Fig. 23. — The same as Fig. @J but for the region of EGOs 03, 05, and 10. 




Fig. 24.— The same as Fig. HJ but for the region of EGOs 07, 09, 11, 13-14, 17-18. The 
location of the Class sources VLA1623-243 is marked with a white pentagram. 



~i — i — I — i — i — i — i — i — I — i — i — i — i — i — I — i — i — i — i — i — I — i — i — i — i — i — I — i — i — r 



22' - 



23 



24' - 



25' - 



-24°26 l 00" I 



■ 



G553Q 




13 
« 



# 11 



14b 




9 



j i i i i i i i i i i i i i i i i i i i i i i i i i i i i i_ 



1 6 h 26 m 27 £ 24 s 



2T 18 a 
RA 



15 E 



Fig. 25. — The same as Fig. [3J but for the VLA1623-243 northwest region. The location 
of the Class sources VLA1623-243 is marked with a pentagram. The GSS source is from 
Grasdalen et al.l (Il973l ). 




Fig. 26 . — The same as Fig. [3 but for the VLA1623-243 southeast region. The OPH source 
is from iGrasdalen et al.l (119731 ) . 
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Fig. 27. — The same as Fig. [3l but for the E GOs 19, 21, and 25 re gion. Source design ations 
used a re as follows: BK LT sources are from iBarsonv et al.l (119971 ). Elias sources fromlElias 
(119781 ) . SR sources from lStruve fc Rudkiobind (119491 ) . VSSG sources from lVrba et al.l (119751 ) 
and WLY sources from IWilking et al. Jl989h . The bow shock morphology of EGO 19 is 
outlined with a black line. 




Fig. 28. — The same as Fig. HJ but for the EGOs 19, 21, and 25 region. 
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Fig. 29. — The same as Fig. El but for the E GOs 20 and 22- 2 3 reg ion. Source designations 
used are as follows: B KLT sources are from IBarsonv et al.l (119971). IS O-oph sources from 
Bontemps et al.l (120011 ) and BBRCG sources from IBarsonv et al.l (119891 ). 
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Fig. 30.— The same as Fig. H but for the EGOs 20 and 22-23 region. 
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Fig. 31. — The sam e as Fig. [3l but for the EGO 24 region. BBRCG sources are from 
Barsony et al.l (119891 ). The bow shock morphology of EGO 24a is outlined with a white line. 




Fig. 32. — The same as Fig. HJ but for the EGO 24 region. 
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Fig. 33. — The sam e as Fig. [31 but for the EGO 26 region. The BKLT source is from 
Barsony et al.l (Il997l ). 




Fig. 35. — The sa me as Fig. [3j but for th e EGOs 27-29 region . YLW sources are from 
Young et al.l (119861 ) and WLY1 sources from IWilking et al.l (Il989l ). 
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Fig. 36. — The same as Fig. HI but for the EGOs 27-29 region. 
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Fig. 37. — The same as Fig. [31 but for the EGOs 28 and 31 region. BKLT sources are from 
Barsonv et al.l (19971). ISO-oph so urces from iBontemps et al.l (120011 ). and AOC sources from 
Alves de Oliveira fc Casalil (b()08h . 




Fig. 38. — The same as Fig. HJ but for the EGOs 28 and 31 region. 
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Fig. 39. — The sam e as Fig. [3j but for the EGOs 30 and 32 region. The Y LW source is from 
Young et al.l (I1986I ) and the GY92 source is from iGreene fc Young! (I1992I ). 
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Fig. 40. — The same as Fig. HI but for the EGOs 30 and 32 region. 
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Fig. 41. — The distribution of IRAC mid-infrared outflows, HH objects, and CO outflows 
in the p Ophiuchi molecular cloud. The grey scale image is the IRAC 4.5 /im image. The 
contours represent antenna te mperatures T* A of 2,4,5,6,7,8,10,12,14,18 and 20K of the 13 CO 
emission. CO outflows are from lWu et al.l (12004 ) and HH objects are from Simbad. The arrow 
marks the direction of the Upper Scorpius OB association with respect to the Ophiuchus 
cloud. The distance shown above the arrow is from the center of current figur e to th e center 
of Upper-Scorpius OB association region that was defined by Ide Zeeuw et al.l (119991 ) . 




Fig. 42. — The distribution of IRAC mid- infrared outflows compared with that of YSO 
candidates. Others are the same as in Fig. SU 




Fig. 43. — The distribution of IRAC mid - infrar ed outflows compared with that of MMSs 
and SMMs. MMSs are from lYoung et al.l (120061 ) and SMMs are from Simbad. Others are 
the same as in Fig. SU 



